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Abstract 

Recently developed time-independent bound-state perturbation theory 
is extended to treat the scattering domain.The changes in the partial wave 
phase shifts are derived explicitly and the results are compared with those 
of other methods. 
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1 Introduction 

In our previous papers [1-6] published recently, a time-independent novel per¬ 
turbation theory has been developed in the bound state domain, which is non- 
perturbative, self-consistent and systematically improvable, and used to treat 
successfully significant problems in different fields of physics. Gaining confidence 
from these applications, we aim through the present work to show that similar 
techniques can also be used in the continuum. 

In the next section we summarize the main ideas of our approach. The ex¬ 
tension of the model for scattering states and the relationship to some other 
perturbation approaches are discussed in Section 3. The paper ends with a brief 
summary and concluding remarks. 

2 The Model 

Let us start with a brief introduction of the formalism to remind the compact 
form of the method, which would provide an easy access of the scheme to under¬ 
standing of the treatments in the continuum. For the consideration of spherically 
symmetric potentials, the corresponding Schrodinger equation in the bound state 
domain for the radial wave function has the form (h = 2m = 1) 
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where Vo is an exactly solvable unperturbed potential together with the angular 
momentum barrier while AV is a perturbing potential. Expressing the wave 
function \P n as a product 

M r ) = X„(r)0 n (r), (2) 

in which \n is the known normalized eigenfunction of the unperturbed Schrodinger 
equation whereas (j> n (r) is a moderating function corresponding to the perturbing 
potential. Substituting (2) into (1) yields 


v" 
A r 


+ — + 2 


Xn <P'n 


Xn 4*n Xn 


= V - E n 


(3) 


Instead of setting the functions \n and 0 n (r), we will set their logarithmic deriva¬ 
tives 


W n = AW n = -fi 
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which leads to 
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where £ n is the eigenvalue of the exactly solvable unperturbed potential, and 


A Wl - A W' n + 2W n AW n = AV(r) - Ae n (6) 

\ rVi Xn 0n / 

in which Ae n is the energy value for the perturbed potential leading to E n = 
e n + Ae n . If the whole potential, involving the perturbing piece AV, can be 
analytically solvable, then Eq.(l) through (5) and (6) reduces to 


(W n + AIE n ) 2 - (W n + A W n y = V - E n , (7) 

which is known as the usual supersymmetric quantum mechanical treatment [7] 
in the literature. 

However, if the whole potential has no analytical solution as the case consid¬ 
ered in this Letter, which means Eq.(6) cannot be exactly solvable for AW, then 
one can expand the functions in terms of the perturbation parameter A, 


AV(r; A) = £ X N AV N (r), AW n (r; A) = £ X N AW nN (r), Ae n (X) = ^ X N Ae nN 

N =1 N= 1 N=1 

( 8 ) 
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where N denotes the perturbation order. Substitution of the above expansion 
into Eq.(6) and equating terms with the same power of A on both sides yields up 
to for instance 0(A 3 ) 

2W n AW nl - A W' nl = m - Ae nl , (9) 

AW* + 2W n AW n2 - A W' n2 = AV 2 - Ae n2 (10) 

2(W n AW n3 + AW nl AW n2 ) - AW' 3 = AK 3 - Ae n3 (11) 

Eq.(6)and its expansion through Eqs.(9-ll) give a flexibility for the easy calcu¬ 
lations of the perturbative corrections to energy and wave functions for the nth 
state of interest through an appropriately chosen perturbed superpotential, ft 
has been shown [1-6] that this feature of the present model leads to a simple 
framework in obtaining the corrections to all states without using complicated 
mathematical pro cedures. 

3 Application to the scattering domain 

It is well known that there are many scattering problems in which the interaction 
between the projectile and the target decomposes naturally into two parts (V = 
Vo + AV). This division is especially useful if the scattering wave function under 
the action one part can be obtained exactly (Vo), while the effect of the other 
(AK) can be treated in some approximation as in the present formalism. 

For simplicity, we here confine ourselves to s —wave scattering from a potential 
which is assumed that vanishes beyond a finite radius R. The associated total 
wavefunction behaves at large distances 

T(r) = — sin (kr + 5), r > R, (12) 

k 

where 5 is the s —wave phase shift. 

Our present treatment of scattering has concerned itself primarily with de¬ 
termining how the solutions of the free Schrodinger equation are affected by the 
presence of the interaction. Within the framework of the present formalism we 
suppose that the solutions of Eq.(5) are known, or are easily found, to give the 
corresponding phase shift 8 0 - Considering the expansion 8 — 5 0 + Ahi + A 2 <5 2 + ..., 
as in Eq.(8), we aim here to derive explicitly solvable and easily accessible ex¬ 
pressions for the phase shift contributions at successive perturbation orders. 

3.1 First-order phase shift correction 

Keeping in mind Eq.(12) and considering the discussion in Section 2, at the first 
perturbation order one has 

(W + AAkVi) = —k cot (kr + <5 0 + A<h), W n = — — = —k cot (kr + <5 0 ), (13) 

X 
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from where the superpotential relating to the perturbing interaction 


AWi(r) 


kSi 

sin 2 (kr + ho) 


( 14 ) 


is obtained assuming that sin XS i = AJi and cosAJi = 1. 

In the second step, one needs to employ Eq. (9) to arrive at another expression 
for AHq. Rearranging the terms, A W[ — 2WAW\ = (Aeq — AVi) and multiply 
both sides by the integrating factor exp(—2 / 0 ' W(z)dz), which is the square of 
the unperturbed wave function x 2 (r) through Eq.(4), one obtain 


= X 2 (r)( A^-AEi). (15) 

The integration, and the remove of Aeq term due to the consideration of elastic 
scattering process here, yields 


U, 

dr* L' 


X 2 (r)AW 1 (r) 


AWi(r) 


~4n / X 2 (-)AEi (z)dz. 
X in Jo 


(16) 


As x — i sin (kr + Jo) i n the asymptotic region, comparison of Eqs.(14) and (16) 
reproduces the first-order change in the phase shift 

r oo 

5i — —k y 2 (r)A Vi(r)dr. (17) 

Jo 

If necessary, the corresponding change in the wavefunction can easily be obtained 
by the substitution of Eq.(16) into (4), </q = exp(— / AW\). For the reliability of 
the present expression obtained, Eq (17), one may compare it with that repro¬ 
duced by other methods. For example, in the limiting case where the unperturbed 
potential vanishes, the unperturbed s— wave function is reduced to a plane wave 
y(r) = sin (kr)/k, and the first-order change in the phase shift becomes 

1 r°° 

Ji = — — / sin 2 (kr)AVi(r) dr (18) 

k Jo 

which is just the first Born approximation for the phase shift [8]. In addition, 
the well known expression for s— wave scattering amplitude by the two-potential 
formula in scattering theory [8], 

roo 

fi = —e 2lSo / x 2 (r)AVi(r)dr (19) 

Jo 

where the phase factor in front of the integration arises because of the standing 
wave boundary conditions, justifies once more our result since fi = —e 2lSo 5i/k 
and, equating this to the above equation leads immediately to Eq.(17). 

The present result has a widespread applicability, which may also be used in 
the treatment of scattering length problems. At low-energy limit, the phase shift 
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is related to the scattering length h^o ~ka where a = ao + Aai + A 2 a 2 + ... 
may be expanded in a perturbation series similar to the phase shift. Outside the 
range of the potential, the unperturbed wave function behaves as x ~* ( r ~ a o )• 
Thus, the hrst correction to the scattering length is 


a i 



ao) 2 AVi(z)dz 


( 20 ) 


which can be calculated for a given AV\. The scattering length has an important 
physical significance. In the low-energy limit only the s —wave makes a nonzero 
contribution to the cross section, so that the angular distribution of the scattering 
is spherically symmetric and the total cross section is 47r(a 0 + Aai + ••O 2 - This is 
also exactly the result obtained in most textbooks for the low-energy scattering 
of a hard sphere of radius Thus the scattering length is the effective radius of the 
target at zero energy. 

As a last example, consider the case of the angular momentum barrier as the 
unperturbed potential Vo = £(£ + 1 )/r 2 that produces [rje(kr)\ with a phase shift 
So = —fbr/2. For a trivial perturbation let us choose AVi = A/r 2 , due to which 
the angular momentum is slightly perturbed I ~ l + A/ (2£ +1) + O (A 2 ). Therefore 
the phase shift correction at first-order is hi = —tt/2(2£ + 1). Again, this exact 
result confirms the reliability of Eq.(17). 


3.2 Second-order phase shift correction 

To solve Eq.(10) for A IT 2 we mimic the preceding calculation. The integration 
factor is the same. In fact, examining Eqs.(9) and (10), the only difference is 
that the quantity AVi — Aeq is replaced by AV 2 — AW 2 — Ae 2 . As Ae 2 term is 
zero due to the process of interest, AIT 2 is thus 


AW 2 (r) 



AW 2 (z) 


A V 2 (z) 


dz. 


( 21 ) 


Bearing in mind that X = \ sin (hr + ho) f° r the region r > R , the second-order 
expansion in the superpotential similar to Eq.(13) provides another expression 
for AIT 2 which is 


AVT 2 (r) 


, 2 cot (kr + h 0 ) 
1 sin 2 (kr + h 0 ) 


+ 


kS 2 

sin 2 (kr + ho) 


( 22 ) 


Comparison of Eqs.(21) and (22), together with the substitution of (14) in (21), 
leads to an auxiliary function for the second order phase shift correction, 


h 2 (r) = [ A V 2 (z) sin 2 (kz + S 0 )dz + /ch 2 [ 

k Jo Jo 


dz 


sin 2 (kz + h 0 ) 


— h 2 cot (kr + ho), 
(23) 
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where a singularity appears in the second integral at z = 0. This problem can 
be circumvented by replacing the lower limit of the integral with R. Assuming 
AV = AVi as in realistic problems of nuclear physics, which means that AV 2 = 0, 
the r —dependent phase shift correction in the second-order is given in the form 
of 

82 (t ) = 8 \ cot (kR + <5 0 ) — 2h 2 cot (kr + 8 0 ). (24) 

As an alternative treatment, which leads to a concrete comparison, one can go 
back to Eq.(21) and split x 2 Ahh 1 2 term in two parts as (x 2 AHA) (AHA) allowing 
to invoke Eq.(16). In this case the comparison of the result with the expansion 
in (22) gives 


5 2 



J z x 2 (y)^yi(y)dy 


81 

k 


cot (kR + 5q) 


(25) 

which is in agreement with the work in [10]. In addition, the use of (17) in 
(24) transforms it into Eq. (25). Furthermore, the reader is reminded that the 
second Born approximation for the phase shift can be most easily derived using 
the variable phase equation approach [11], 


poo pr 

82 = 2 k 2 / x 2 (r)AEi(r) cot (kr)dr / ~x 2 (y)AVi(y)dy 
J 0 Jo 


(26) 


which, in the light of Eq. (15), is the same result as we find from Eq (25), by 
putting <5 0 = 0 . Higher order terms can also be evaluated in the same manner. 


4 Concluding Remarks 

The recently introduced time-independent perturbation theory has been success¬ 
fully extended from the bound state region to the scattering domain. For the 
clarification, the work has been carried out with the consideration of s —wave 
scattering only. However, generalization of the formalism to higher partial waves 
in the scattering domain does not cause any problem. The inclusion of the cen¬ 
trifugal barrier contribution in the effective potential for instance leads to the 
replacement of the s —wave phase shift with Si — £n /2 due to the related wave 
function y(r) = sin (At + 81 — iTr/ 2 )/k in the asymptotic region, supposing both 
the unperturbed and perturbed potentials vanish at a large r > Ri which means 
that in the region R\ < r < R there is then only the centrifugal barrier contri¬ 
bution. This inclusion requires simply to repeat the present calculations for the 
replacement in the phase shift. 

It should be stress that, anything that can be achieved from the present 
formalism must also be obtainable from the works [9,10] in the literature. For 
instance, considering the bound state region, Bender’s formalism [9] can be simpli¬ 
fied by introducing the auxiliary function i7v(r) such that the whole wave function 
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Tjv(r) = x(r)lbv(r) where denotes the perturbation order. The first-order cor¬ 
rection can then be written as 4- X 2 ^£r — (AVi — AeQx 2 which corresponds ex¬ 
actly to the present treatment by Eq. (15) when we identify AWi = dF/dr. The 
higher order calculations can be linked to ours in the similar manner. Whereas, 
the works of Milward and Wilkin [10] may be related to the present formalism 
in both domain, the bound and scattering region by making a relation between 
their probability density distributions/derivatives and our AW functions, such 
as AWq = —Pq/2P 0 at the zeroth order, AWi = {—Pi/^Po)' at the first order 
and Al'l'b = (—P2/2P o y at the second order etc. Nevertheless, the present tech¬ 
nique provides a clean and explicit route for the calculations without tedious and 
cumbersome integrals. 

The energy variation of the scattering wave function and phase shift can also 
be studied by perturbing in the energy. We wish to stress that all these effects 
depend purely upon the perturbation and the unperturbed wave function; explicit 
knowledge of the unperturbed potential is not necessary. This exposition will be 
deferred to a later publication. 
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